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I.  INTRODUCTION 


Because  of  failures  in  structures  caused  by  ground  motion  from  earth¬ 
quakes  or  explosives,  soi I -concrete  interaction  problems  are  playing  an 
increasingly  important  part  in  engineering  design  procedures  and  safety  analy¬ 
ses.  The  information  available,  however,  is  still  sparse  due  to  the  complex¬ 
ity  of  the  nonlinear  response  of  these  geologic  materials. 


In  the  past,  many  models  for  representing  interfaces  were  based  on 
Coulomb's  law,  i.e.,  no  slip  between  two  bodies  will  occur  unless  a  critical 
tangential  stress  is  reached,  and  the  relative  motion  is  constrained  to  be 
along  the  contact  surface  if  the  critical  state  is  reached.  For  dissimilar 
materials  like  concrete  and  soil,  however,  Coulomb's  law  may  not  be  appli¬ 
cable,  since  slip  does  not  occur  at  the  interface.  Instead  of  slip,  a  band 
with  large  shear  deformation  is  observed  in  the  weaker  of  the  two  materials, 
and  these  observed  zones  with  large  strain  gradients  cannot  be  automatically 
predicted.  Nevertheless,  the  simple  friction  approach  has  been  utilized  with 
either  slip  lines  or  special  interface  elements  to  obtain  computational  pro¬ 
cedures,  which  in  turn  introduce  questions  of  stability  and  computational 
complexity. 


This  research  addresses  a  new  approach  for  simulating  a  soi 1 -concrete 
interface  where  both  theoretical  and  computational  aspects  are  kept  in  mind. 
The  idea  relies  on  the  observations  that  (1)  slip  does  not  occur  for  the  soil  - 
concrete  interface  until  a  failure  state  is  reached,  and  that  (2)  progressive, 
distributed  damage  produces  a  band  of  permanent  deformation  in  the  soil  adja¬ 
cent  to  the  interface.  The  lateral  dimension  of  the  band  appears  to  depend  on 
the  physical  characteristics  of  the  materials.  The  only  way  to  predict  a 
finite  zone  of  localized  deformation  under  a  homogeneous  state  of  stress  is  to 
utilize  a  nonlocal  constitutive  law.  Furthermore,  it  appears  that  localiza¬ 
tion  occurs  only  when  softening  is  present.  Thus,  if  a  nonlocal  constitutive 
model  is  postulated  where  softening  is  governed  by  both  strain  and  the  gradi¬ 
ent  of  strain,  then  softening  localizes  into  a  finite  region,  which  cannot  be 
predicted  by  using  a  local  constitutive  model.  With  the  use  of  a  nonlocal 
constitutive  model  for  soil,  and  the  assumption  that  the  interface  initiates 
softening  in  the  weaker  of  the  two  materials,  the  deformation  field  adjacent 
to  the  interface  can  be  simulated. 


In  this  study,  numerical  solutions  of  a  model  problem  with  dynamic  and 
static  shear  loadings  are  compared  with  other  numerical  data.-  It  was  conclude 
that  the  new  approach  is  reasonable  in  an  engineering  sense,  and  it  is  feasi¬ 
ble  to  incorporate  the  model  into  existing  finite  element  and  finite  differ* 

ence  computer  codes. 


II,  SCOPE  OF  THE  INVESTIGATION 


OBJECTIVE 

Based  on  existing  experimental  data,  the  general  behavior  of  soil  adja¬ 
cent  to  a  concrete  interface  can  be  described.  Some  assumptions  are  proposed 
to  construct  a  new  nonlocal  constitutive  model  for  simulating  the  behavior  of 
soil  in  the  vicinity  of  an  interface.  Using  computational  plasticity  theory, 
the  characteristics  of  the  model  are  discussed.  Existing  theoretical  data 
are  used  for  comparison  with  the  model.  Conclusions  and  conments  on  possible 
future  research  are  given, 

INTERFACE  PHENOMENA 

Interface  behavior  has  drawn  considerable  attention  in  recent  years, 
especially  for  the  following  engineering  problems:  retaining  wall  design, 
skin  friction  in  piles,  soi 1 -structure  interaction  for  nuclear  power  plants 
subjected  to  earthquake  loadings,  and  for  underground  structures  subjected  to 
blast  and  ground  shock. 

Due  to  the  limitations  of  experimental  techniques,  and  the  assumptions 
that  Coulomb's  friction  principle  holds  for  interfaces  of  deformable  bodies, 
there  have  been  a  limited  number  of  investigations  of  interface  phenomena. 
Existing  experimental  data  do  not  provide  an  accurate  view  of  the  deformation 
field  in  the  vicinity  of  the  interface.  Nevertheless,  an  approximate  descrip¬ 
tion  and  some  knowledge  of  the  soi 1 -concrete  interface  can  be  obtained  by 
carefully  analyzing  the  data.  This  research  discusses  the  problem  from  a 
continuum  point  of  view  and  presents  an  alternate  viewpoint  based  on  the 
observation  that,  adjacent  to  an  interface,  a  region  with  a  large  gradient  of 
strain  is  often  formed  for  shear  loading  (Fig.  1).  An  important  part  of  this 
work  was  to  develop  a  procedure  for  describing  this  phenomenon. 

Among  the  investigators  dealing  with  soi  1 -concrete  interface  phenomena, 
Huck,  et  a  1  - ,  (Refs.  1  and  2),  were  unique  in  that  both  experimental  and  ana¬ 
lytical  investigations  were  performed.  Also,  several  valuable  references  were 
given.  From  their  interface  tests  under  static  and  dynamic  shear  loadings,  it 
was  concluded  that  interface  shear  resistance  is  strongly  influenced  by  the 
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Figure  1.  Geometry  for  the  interface  problem. 
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soil  type,  the  stress  normal  to  the  interface,  and  the  relative  size  of  the 
interface  roughness  with  respect  to  the  soil  grain  diameter.  Total  interface 
shear  resistance  is  reduced  below  the  soil  shear  strength  due  to  softening  of 
the  expansive  soil  and  loss  of  adhesion.  It  appears  that  localization  of 
large  deformations  usually  occurs  adjacent  to  the  interface,  in  a  region  char¬ 
acterized  as  a  shear  band,  which  increases  gradually  in  thickness  with  contin¬ 
ued  increase  in  applied  static  displacement.  After  each  interface  test,  the 
final  dimension  of  the  band  of  permanent  deformation  Is  not  the  same  for 
different  materials. 

Drumm  and  Desai  (Ref.  3)  worked  on  sol  1 -structure  interface  tests  under 
static  and  cyclic  loadings  by  using  a  now  testing  apparatus  called  the  "Cyclic 
Multi -Degree-of-Freedom"  device,  They  indicated  that  the  behavior  of  a  struc¬ 
ture  resting  on,  or  embedded  in  a  foundation  depends  on  the  material  proper¬ 
ties  of  both  the  structure  and  the  foundation  medium,  and  on  the  types  of 
applied  loadings.  Soi 1 -structure  Interaction  becomes  more  significant  when 
structures  are  subjected  to  dynamic  loadings.  The  interface  strength  was 
found  to  be  less  than  that  of  soil,  and  gradual  separation  may  occur  at  or 
near  the  interface, 

Vermeer  (Ref.  4)  investigated  interface  behavior  between  granular  materi¬ 
als  and  metals  during  powder  compaction  in  a  die.  There  are  two  friction 
angles  in  powder  compaction:  one  for  powder-wall  friction,  and  another  for 
powder-powder  friction.  These  two  angles  were  shown  to  be  related  for  compac¬ 
tion  in  a  die  with  rough  walls.  From  tests,  it  was  pointed  out  that  the 
powder-wall  slip  plane  is  In  fact  a  thin  shear  band. 

A  number  of  factors,  not  directly  pertinent  to  the  basic  soi 1 -concrete 
interface  problem,  might  act  to  obscure  the  true  phenomena.  However,  by  ana¬ 
lyzing  the  experimental  apparatus  and  test  results  mentioned  above,  it  was 
possible  to  draw  the  following  descriptions  for  the  mechanism  of  an  interface 
consisting  of  concrete  and  particulate  material. 

Imperfections  at  or  near  the  soi 1 -concrete  interface  will  initiate  pro¬ 
gressive  distributed  damage  such  as  dispersed  microcracking,  void  formation, 
or  loss  of  interparticle  contacts.  From  the  macroscopic  viewpoint,  large 
deformations  can  be  observed  adjacent  to  the  interface,  and  the  size  of 


the  band  of  permanent  deformation  will  enlarge  with  Increase  in  the  applied 
displacement.  Although  relative  motion  occurs  between  particles  as  soon  as 
external  loading  is  applied,  no  slip  seems  to  occur  as  the  stress  is 
increased.  If  deformation  continues  to  increase,  the  stress  will  reach  a  peak 
and  then  decrease  in  what  Is  called  the  postpeak  regime.  Slip  appears  to  be 
initiated  somewhere  in  the  postpeak  regime.  When  macroscopic  slip  occurs,  the 
dimension  of  the  band  seems  fixed  and  the  total  deformation  field  is  not  con¬ 
tinuous,  Thus,  slip  might  be  considered  as  a  mechanism  that  occurs  after  the 
formation  of  the  shear  band. 

PRESENT  STATUS  OF  MODELS  FOR  INTERFACES 

During  the  past  20  years,  evaluation  of  the  dynamic  interaction  between 
a  structure  and  the  surrounding  medium  has  been  accomplished  by  two 
approaches,  consisting  of  either  the  use  of  classical  continuum  mechanics,  or 
the  finite  element  method. 

The  first  approach  simulates  the  effects  of  soil  on  the  structural 
response  by  a  series  of  springs  and  dashpots  representing  a  theoretical  half 
space  surrounding  the  structure.  This  method  is  generally  limited  to  elastic 
or  viscoelastic  representations  of  soils.  The  finite  element  method  is  an 
alternative  approach  to  solving  problems  In  continuum  mechanics.  Although 
based  on  the  same  principles  of  mechanics,  it  changes  the  partial  differential 
equations  governing  the  motion  of  the  continuum  Into  a  finite  set  of  ordinary 
differential  equations  in  the  time  domain  which  makes  it  possible  to  incorpo¬ 
rate  a  more  reasonable  constitutive  relation  into  the  interface  model.  A 
survey  of  both  the  spring-dashpot  and  finite  element  approaches  to  the 
structure-soil  interaction  problem  was  given  by  Madjic\n,  et  al.,  (Ref.  5). 
Also,  Seed,  et  al.,  (Ref.  6),  presented  an  evaluation  of  the  strengths  and 
weaknesses  of  the  two  approaches.  Although  the  above  information  reviewed 
only  the  efforts  In  an  infant  period  of  research,  the  history  of  development 
in  this  field  was  quite  helpful. 


Recently,  emphasis  in  research  nas  been  on  improving  the  description  of 
interface  behavior,  not.  only  in  the  prepeak  regime,  but  also  in  the  postpeak 
regime,  while  computer-based  techniques  have  supplied  powerful  tools  for 
analyses.  Mathematical  models  are  set  up  witn  the  use  of  both  conservation 
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principles  and  constitutive  relations.  Conservation  principles,  based  upon 
classical  physics,  or  Newtonian  physics,  have  been  well  established.  Hence, 
accurate  descriptions  of  the  behavior  of  materials  have  become  more  and  more 
important.  Desai  and  Siriwardane  (Ref.  7)  reviewed  and  explored  the  constitu¬ 
tive  behaviors  of  engineering  materials,  with  emphasis  on  geolo9-cal  materials 
such  as  soil  and  concrete.  Currently,  in  addition  to  the  classical  continuum 
approach  for  soi 1 -concrete  interface,  both  the  particle  approach,  together 
with  statistical  analysis,  and  the  nonlocal  continuum  approach  have  been 
proposed. 

Huck,  et  al . ,  (Refs.  1  and  2),  derived  a  soi 1 -concrete  interface  model 
using  basic  frictional  mechanisms  as  applied  to  a  particulate  soil.  The  par¬ 
ticle  approach  employed  the  mechanisms  of  adhesion,  plowing,  and  lifting.  The 
concrete  interface  surface  was  considered  as  a  flat  plane,  with  protruding 
spherical  caps  representing  the  roughness  or  asperities  cf  the  concrete  sur¬ 
face.  Soil  in  contact  with  the  concrete  interface  was  modeled  as  spherical 
particles  of  known  size,  but  randomly  distributed  over  the  interface  surface. 
In  the  model,  relative  slip  between  soil  and  concrete  occurred  if,  and  only 
if,  a  critical  tangential  stress  was  reached,  and  the  relative  motion  was 
assumed  to  be  along  the  contact  surface.  The  tangential  stress  was  propor¬ 
tional  to  the  normal  stress.  It  appears  that  gradual  formation  of  a  shear 
band  could  not  be  predicted.  This  constitutive  model  is,  by  nature,  a  classi¬ 
cal  Mohr-Coulomb  model. 

Vermeer  (Ref.  4)  investigated  interface  response  in  powder  ;ompaction  by 
means  of  a  nonassociated  perfect-plasticity  model.  The  idea  of  sliding  on  a 
single  plane  was  considered  much  too  restrictive.  Therefore,  a  new  idea 
allowed  small  slipping  motions  on  an  infinite  number  of  closely  spaced  paral¬ 
lel  surfaces,  and  used  the  Mohr-Coulomb  criterion.  To  calculate  the  stress 
and  density  distributions  at  several  stages  of  the  loading  process,  three 
approaches  were  suggested.  The  first  was  the  use  of  a  nonassociated  consti¬ 
tutive  model  and  a  fine  mesh  near  the  walls  of  the  die,  so  that  a  thin  shear 
band  could  develop.  The  second  consisted  of  special  thin  interface  elements 
to  model  the  shear  band  along  the  wall.  Frictional  sliding  was  implemented 
directly  as  a  mixed  boundary  condition  in  the  third  approach. 
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Drumm  and  Desai  (Ref.  3),  and  Oden  and  Pi  res  (Ref.  8),  describe  the  con¬ 
stitutive  behavior  of  the  interface  by  using  frictional  principles  with  some 
modifications.  Either  elasticity  or  perfectly-plastic  relations  have  been 
used  for  simulating  soi 1 -concrete  interface  by  these  authors  as  wel I  as 
others . 

From  the  research  efforts  mentioned  above,  it  was  found  that  assumptions 
based  upon  frictional  theories  were  commonly  adopted,  and  numerical  methods 
are  under  active  development  to  solve  the  continuum  models.  It  is  doubtful, 
however,  that  frictional  theories  are  useful  for  representing  interfaces 
involving  materials  with  disparate  properites  such  as  concrete  and  soil.  In 
addition,  there  are  some  unpleasant  features  associated  with  numerical  proce¬ 
dures  for  frictional  models,  such  as  size-dependent  response  and  computational 
difficulties  due  to  special  algorithms.  To  obtain  a  basic  understanding  of 
the  inherent  physical  nature  of  the  interface,  and  to  apply  the  knowledge  to 
engineering  practice,  the  exploration  of  a  more  reasonable  and  convenient 
approach  seems  necessary. 

With  the  development  of  knowledge  for  heterogeneous  materials,  a  nonlocal 
continuum  theory  offers  hope  for  modeling  the  postpeak  behavior  of  geological 
materials.  The  word  "nonlocal"  in  this  context  means  that  stress  depends  both 
on  strain  and  strain  distribution.  The  latter  can  be  given  through  the  use  of 
integrals  and  a  weighting  function,  or  by  derivatives  of  strain.  The  full 
force-deformation  behavior,  but  not  necessarily  the  stress-strain  relation¬ 
ship,  beyond  peak  stress  has  been  documented  for  concrete,  rock,  soils,  and 
sor.ie  other  materials  under  uniaxial  stress.  One  approach  was  to  model  the 
postpeak  behavior  by  using  plasticity  and  the  reduction  of  stress  with  strain. 
However,  serious  theoretical  and  computational  difficulties  caused  by  strain 
softening  have  occurred,  and  they  have  been  intensely  debated  at  recent  con¬ 
ferences.  Valanis  (Ref.  9),  Sandler  and  Wright  (Ref.  10),  and  Read  and 
Hegemier  (Ref.  11)  discussed  in  detail  the  various  arguments  concerning  strain 
softening.  Bazant,  et  al . ,  (Ref.  12),  incorporated  strain  softening  into  a 
consitutive  model  by  using  a  nonlocal  continuum  theory.  Schreyer  and  Chen 
(Ref.  13)  investigated  the  effects  of  strain  softening  and  localization  in  a 
structure  by  assuming  that  stress  was  a  function  of  both  strain  and  strain 
gradients.  It  was  concluded  that  current  experimental  techniques  could  not 
rule  out  strain  softening  because  instrume  itation  has  not  been  designed  to 
identify  localization  and  to  measure  strains  in  the  resulting  small  zones. 


8 


Also,  several  types  of  nonlocal  approaches  have  overcome  obstacles  caused  by 
strain  softening.  At  present,  (1)  a  finite  strain-softening  zone,  (2)  nonzero 
energy  dissipation,  and  (3)  numerical  stability  of  dynamic  problems  have  been 
achieved.  Additional  research  is  required  to  generalize  the  theory  on  strain 
softening.  Incorporating  strain  softening  into  constitutive  models  would  be 
useful  and  feasible  in  engineering  practice  if  the  response  of  a  material  can 
be  predicted. 

Considering  the  interface  phenomena,  it  has  been  suggested  that  a  one¬ 
dimensional  model  (Ref.  13)  be  used  to  simulate  the  weaker  of  two  materials  at 
an  interface,  with  the  assumption  that  an  interface,  in  effect,  weakens  the 
softer  material.  If  an  interface  can  be  dealt  with  in  such  a  manner,  then 
special  algorithms  would  not  be  necessary  to  handle  interfaces,  and  further¬ 
more,  a  band  of  permanent  deformation  adjacent  to  the  interface  will  be  pre¬ 
dicted  which  is  representati ve  of  physical  behavior.  Using  conventional , 
constitutive-equation  algorithms  and  finite-element  methods,  static  and 
dynamic  response  features  of  the  interface  can  be  obtained  with  a  simple, 
convenient  calculation. 

Since  an  ultimate  objective  i •  *o  have  a  completely  general  model,  this 
research  provides  a  three-dimen  al  formulation.  However,  applications  are 
limited  to  the  one-dimensiona1  i. 

PROCEDURE 

Several  assumotio  ;osed  based  upon  existing  experimental  data. 

A  three-dimensionat  .  oi.stitutive  model  is  then  developed  for  simu¬ 

lating  the  soil -concrete  interface  by  the  use  of  nonlinear  plasticity  theory. 
The  model  reflects  the  characteristics  of  geological  materials  more  accurately 
than  a  bilinear  constitutive  model.  Features  of  the  model  are  explored  under 
different  values  of  mean  pressure  and  strain  gradient.  In  order  to  verify  the 
model  for  an  interface  under  dynamic  shear  loading,  the  finite  element  method 
is  used,  together  with  explicit  time  integration  and  coinput  tional  plasticity 
techniques,  to  obtain  numerical  solutions  which  are  compared  with  existing 
theoretical  data  for  the  wave  propagation  problem.  The  conclusion  of  this 
research  shows  results  for  static  she.  oading. 
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III.  MATHEMATICAL  MODEL 

ASSUMPTIONS 

Making  assumptions  was  the  first  nontrivial  step  in  the  formulation  of 
the  new  model.  Reasonable  assumptions  yield  good  results  with  a  modest 
effort,  whereas  poor  assumptions  can  lead  to  considerable  computational 
effort.  Making  good  assumptions  depends  on  engineering  experience  and  intui¬ 
tion.  Furthermore,  each  assumption  has  its  own  limitation  which  must  corre¬ 
spond  to  the  problems  being  considered. 

As  mentioned  previously,  there  is  not  enough  experimental  data  to  provide 
a  good  qualitative  description  of  the  behavior  at  soil -concrete  interfaces. 

For  example,  it  is  still  uncertain  at  what  stage  of  the  loading  path  the  mac¬ 
roscopic  slip  between  concrete  and  soil  occurs,  if  it  occurs  at  all.  The 
following  assumptions  were  made  to  construct  a  new  nonlocal  continuum  model, 
which  is  different  from  the  usual  Mohr-Coulomb  model. 

Assumption  1--Slip  does  not  occur  at  the  contact  surface  between  concrete 
and  soil  under  shear  loading.  Instead,  with  an  Increase  in  load,  a  shear  band 
will  form  due  to  dispersed  microcracking.  This  shear  band  represents  a  local¬ 
ization  of  deformation  that  is  accompanied  by  softening. 

Assumption  2— The  imperfection  of  the  soi 1 -concrete  interface  reduces  the 
strength  of  the  soil,  so  that  with  large  enough  loads,  the  softening  is  initi¬ 
ated  in  the  region  of  soil  adjacent  to  the  interface.  With  this  assumption, 
the  roughness  of  the  concrete  surface  is  simulated  through  the  strength  prop¬ 
erty  of  the  adjacent  material,  rather  than  through  a  friction  coefficient. 

The  result  is  that  softening  is  initiated  at  a  lower  load  for  a  smooth  inter¬ 
face  compared  to  a  rough  interface,  but  the  evolution  of  softening  and  the 
localization  zone  will  be  governed  by  the  constitutive  equation. 

Assumption  3- -A  nonlocal  constitutive  theory  is  required  to  represent,  the 
localized  behavior  of  softening  in  soil.  This  theory  provides  a  means  for 
controlling  the  size  of  the  localization. 

Assumption  4--A  suitable  nonlocal  theory  is  obtained  if  the  limit  stress 
depends  on  the  strain  gradient.  Read  and  Hegemier  (Ref.  11)  indicated  that 
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strain  softening  is  not  a  true  material  property,  but  is  simply  a  manifesta¬ 
tion  of  the  effects  of  progressi vely  increasing  inhomogeneity  of  deformation. 
However,  an  alternate  interpretation,  that  equally  explains  the  observed  phe¬ 
nomena,  is  that  softening  occurs  only  in  the  presence  of  strain  gradients. 

The  assumption  that  limit  stress  reduces  with  the  strain  gradient  is  therefore 
a  simple  mathematical  representation  of  the  observation. 

BASIC  EQUATIONS 

To  simulate  the  general  constitutive  behavior  of  soi  1 -concrete  inter¬ 
faces,  a  nonassociated  Prager-Drucker  model,  with  some  modifications,  was 
adopted  to  provide  appropriate  inelastic  response  features  of  soils.  Such  a 
model  does  not  accurately  reflect  the  behavior  of  soils,  but  it  is  one  of  the 
simplest  inelastic  models  that  can  be  used.  If  strain  gradients  are  incorpo¬ 
rated  into  a  complex  model,  it  becomes  difficult  to  separate  and  study  the 
various  features  individually. 

In  the  theory  of  plasticity  (Ref.  14),  the  case  of  combined  stresses  is 
represented  by  means  of  the  concept  of  a  yield  surface.  For  any  point  in  a 
material  region,  assume  that  there  exists  a  domain  in  the  stress  space  in 
which  material  behavior  is  governed  by  elasticity  theory.  The  hypersurface 
bounding  this  domain  is  known  as  the  yield  surface.  If  the  material  is  per¬ 
fectly  plastic,  the  yield  surface  does  not  move.  If  the  yield  surface  expands 
out  to  a  limit  state,  the  material  is  said  to  harden,  while  the  material 
starts  to  soften  as  soon  as  the  yield  surface  contracts  in  a  controlled  man¬ 
ner.  The  motion  of  the  surface  is  governed  by  internal  variables  which  record 
that  part  of  the  stress  history  involving  inelastic  behavior  and  other  varia¬ 
bles  directly  affecting  yield.  The  motion  of  the  stress  point  inside  the 
surface,  i.e.,  in  the  elastic  domain,  is  not  dependent  on  internal  variables. 
In  the  general  constitutive  model  used  in  this  research,  the  internal  vari¬ 
ables  considered  included  an  inelastic  strain  invariant  e1 ,  the  mean  pressure 
P,  and  a  strain  gradient  invariant  g1 .  The  indicial  notation  was  preferred, 
because  this  notation  allowed  us  to  write  equations  in  a  much  shorter  form 
than  would  otherwise  be  possible.  In  representations  of  first  and  second 
order  tensors,  subscripts  refer  to  the  indicial  notation,  while  superscripts 
imply  the  definitions  of  tensors,  e.g.,  the  elastic  or  inelastic  part  of  the 
strain  tensor. 
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In  the  following  part  of  this  subsection,  the  modified  Prager-Drucker 
yield-surface  function  is  discussed  first;  then  the  potential  function  and 
some  important  formulae  are  given.  To  solve  the  nonlinear  constitutive  equa¬ 
tion,  an  incremental  numerical  procedure  is  described  based  upon  computational 
plasticity  theory. 

For  convenience  in  numerical  calculations,  the  yield  surface  is  defined 
in  a  dimensionless  form  as: 


o  -  (1  +  a.P)  G  H 

F  -  - - - —  =  0 

ffLS 


(1) 


in  which  is  the  limit  value  of  the  shear  stress-shear  strain  curve,  and  o 


is  the  second  invariant  of  deviatoric  stress  as  follows: 
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The  function  H  incorporates  both  strain  hardening  and  softening  according  to 
the  following  relations: 


o  < 
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e  <  eL 
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e1  >  (5) 


The  total  strain  tensor  is  divided  into  an  elastic  part  and  an  inelastic  part, 
as  follows: 


Then,  one  of  the  internal  variables  in  Eqs.  4  and  5,  the  inelastic  strain 
i n variant  is  : 


51 
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where 


i  1  1 
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The  values  of  H  at  initial  yield  and  the  limit  state  are  denoted  by  L  and  l 

"•i  •  i  ^ 

respectively.  The  value  of  e  at  the  limit  state,  e^,  is  assumed  to  be  inde¬ 
pendent  of  path.  The  parameter  n  controls  the  shape  of  the  stress-strain 
relation  in  the  hardening  regime  while  a5  controls  the  shape  during  softening 
The  softening  formulation  for  H  given  in  Eq.  5  provides  a  horizontal  slope  at 
the  limit  state  and  a  limiting  value  for  H  of  a4  for  large  e1 . 


The  invariant  (g1 )  of  the  strain  gradient  is  defined  by  the  relation: 


The  strain-gradient  factor,  Gp,  in  Eq.  I  is  assumed  to  be  of  the  form 

Gr  =  (1  -  a^e"^1  +  a2  (10) 

where  the  parameter  a,  controls  the  rate  at  which  G  decreases  with  g1.  For 
— i  — i  ^ 

g  equal  to  zero,  G^  is  unity,  and  for  large  g  ,  the  function  Gr  asymptot¬ 
ically  approaches  a2. 

To  be  consistent  with  experimental  data  (Refs.  1,  3,  and  7)  for 
geological  and  other  frictional  materials,  the  mean  pressure 

P  -  -  -  a,.-  (ID 

3  11 

was  also  incorporated  into  the  constitutive  model  as  a  coefficient  of  the 
parameter  a  1 . 
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With  the  use  of  a  nonassociated  flow  rule,  the  inelastic  part  of  the 
strain  tensor  is  given  by 
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de  •  -  dx  „ — 
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where  the  potential  function  was  chosen  to  be 


(12) 
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(13) 


so  that  the  inelastic  strain  tensor  contains  no  volumetric  component.  This  is 
consistent  with  metal  plasticity  but  is  not  entirely  appropriate  for  soils 
which  display  inelastic  volumetric  strain.  The  inclusion  of  inelastic  volu¬ 
metric  strain  would  have  provided  a  complication  that  has  no  direct  relation¬ 
ship  to  softening  and  localization. 


The  elastic  part  of  the  strain  tensor  is  related  to  the  stress  tensor 
through  the  elastic  constitutive  equation: 


a 


1j 


(14) 


where  the  tensor  is  normally  restricted  to  the  isotropic  case. 

Now  the  general  constitutive  model  is  established.  This  model  is  linear 
in  the  elastic  regime,  and  nonlinear  in  the  hardening  and  softening  regimes. 
Mean  pressure  increases  the  strength  of  the  material,  while  strain  gradient 
reduces  its  strength. 


In  the  numerical  procedure,  an  incremental  scheme  was  adopted  to  solve 
the  nonlinear  constitutive  model.  The  incremental  relations  are  summarized 
below: 
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In  practice,  the  total  strain  increment  Ae. . ,  the  total  strain  from  the 

*  v 

previous  step  (e^.)^,  and  the  stress  from  the  previous  step  (aij)p  are  known. 
The  difficult  problem  is  that  the  elastic  and  Inelastic  part  of  Ae.  .  are 

1  v  j 

unknown.  To  solve  the  above  nonlinear  equations  with  unknowns  e^.,  a^,  e^., 
aa,  and  Aejj,  an  iterative  procedure  was  used.  The  basic  ideas  are  those  of 
computational  plasticity  (Ref.  14).  The  following  is  an  outline  of  the  con¬ 
stitutive  equation  algorithm. 

Step  1— The  user  inputs  components  of  the  total  strain  increment  tensor, 
material  parameters,  and  other  control  parameters.  Existing  values  of  stress, 
total  strain,  elastic  strain  a.id  inelastic  strain  tensors  are  stored  as 

(°i J > P*  (eijV  (eijV  and  (eijV  resPectiye1y- 

Step  2— First,  assume  that  the  step  is  elastic: 


ae .  .  =  Ae .  . 
ij  U 


Step  3— For  future  use,  calculate: 
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I  *  0  (I  =  iteration  count) 


15 


Step  4 — Find  increments  in  stress  components: 

Q 

A<? .  .  a  C .  .  i  i  Ae>  i 

ij  ijkl  kl 

Step  5--Update  old  stress  components  and  the  iteration  count: 


Step  6--Evaluate  o  and  the  yield  function  F. 


Step  7 — Check  the  yield  function  F: 

(a)  If  F  <  0  and  I  =l  (i.e.,  the  first  step),  then  the  solution  is 
elastic  and  the  stresses  are  correct.  Go  to  Step  12. 

Otherwise, 

(b)  If  If  I  <  e  and  I  >  1,  then  the  plasticity  solution  has  been  obtained 
and  the  stresses  are  correct.  Go  to  Step  12. 

Otherwise,  go  to  Step  8. 

Step  8— The  iteration  begins.  The  following  procedure  is  used  to  solve 
for  ax: 


For  the  first  step,  (I  =  1): 


(0 1 

(AX)t  =  -  (one-step  Newton-Raphson  procedure) 

<F,)| 

For  other  steps  (I  >  1): 
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(F)i  '  (F)i-i 


(backward  finite  difference) 


(ax)j  = 


(aa)m  (F)j 
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Step  9--Update  inelastic  strain  increments: 


( Aei  j )  I  "  (Aeij)l-l  *  Tij 

Step  10— Recalculate  elastic  strain  increments: 

Step  ll--Go  to  Step  4. 


Step  12 — Update  inelastic  strain  tensor: 


e1.. 

ij 


Ae 


ij 


Step  13--Compute  inelastic  strain  invariants: 


<-i  ,-i,  -i 

e  =  (e  )p  +  ae 


Step  14 — Update  elastic  strain  and  total  strain  tensors: 


e  ,  e  v  e 

eij  ■  <eij>P  +  4eij 
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e .  .  *  ( e . . )  +  ae .  . 
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Step  15~Exit  for  next  increments  of  total  strain  tensor  components. 

In  the  algorithm  described  above,  the  accuracy  and  number  of  iterations 
are  controlled  through  the  magnitudes  of  the  convergence  parameter  e  and  the 
increments  of  total  strain.  Small  values  of  e  require  generally  more  itera¬ 
tions  and,  hence,  more  computer  time.  This  algorithm  is  simple  and  appears  to 
be  effective  in  an  engineering  sense. 


A  CONSTITUTIVE  RELATION  FOR  THE  INTERFACE  UNDER  SHEAR 


The  constitutive  model  needs  only  to  be  one-dimensional  for  the  soil- 
concrete  interface  under  shear  loading  as  shown  in  Fig.  1,  and  can  be  derived 
directly  from  the  above  general  equations.  The  stress  state  of  the  problem 
is: 


<*11  5  022  *  033  ~  023  s  013 
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The  use  of  Eq.  19  results  in  the  relations: 
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To  simulate  the  behavior  of  soil-concrete  interfaces,  Eqs.  20-23  are  used 
to  represent  the  constitutive  characteristics  of  soil,  and  the  interface  is 
assumed  to  initiate  failure.  The  following  subsection  illustrates  the  effects 
of  mean  pressure  and  strain  gradients  on  the  model  through  the  use  of  the 
constitutive  equation  subroutine.  The  results  should  be  helpful  for 
understanding  the  global  behavior  of  a  soi 1 -concrete  interface. 


FEATURES  OF  THE  MODEL 

A  total  of  eight  parameters  in  Eq.  22  can  be  used  to  simulate  different 

types  of  materials.  Among  these  parameters,  lo  and  Ip  are  the  elastic  limit 
and  the  peak  value  of  the  shear  stress-strain  curve  respectively,  while  a,  and 
a  are  chosen  so  stress  will  not  be  null  in  the  postpeak  regime,  a  reasonable 
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assumption  based  on  experimental  observations.  Parameters  a^  a,,  a?,  and  n, 
control  the  basic  shape  of  the  stress-strain  curve.  For  the  purpose  of 
illustration,  the  effects  of  these  four  parameters  are  shown  separately  in 
Figures  2,  3,  4,  and  5  in  which  stress  has  been  normalized  with  respect  to  the 
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Shear  strain. 


Figure  5.  Effect  of  the  parameter  a*,  or  the  invariant  of  strain 
gradient,  on  stress-strain  response* 


stress  at  the  limit  state,  aLS,  and  strain  has  been  normalized  with  respect  to 
the  strain  at  the  elastic  limit,  e^.  The  other  material  and  control  parame¬ 
ters  used  in  the  numerical  investigations  were  assigned  the  values: 

L0  =  50  /I  Lp  =  100  PT  a2  =  0.1  aL$  =1  g1  =  10 

a,  *  10  fT  - - -  =  2  e\  *  10  ecl  *  25  P  *  100 

2(1  +  v)  L  EL 

e  =  10-**  aLS  ■  100  fZ~  Ae12  =  0.5 

Figure  2  shows  the  effect  of  the  parameter  n  in  the  hardening  regime  for 
ax  «  0,  a3  =  0,  and  a5  =  1.  Note  that  the  curve  is  smoother  with  a  smaller 
value  for  n. 

Figure  3  shows  the  effect  of  the  parameter  as  in  the  softening  regime 
for  a3  =  0,  a3  =  0,  and  n  =  0.5.  By  changing  the  values  of  a5,  different 
postpeak  behaviors  can  be  simulated.  As  a  limiting  case,  the  material  is 
perfectly  plastic  when  a§  is  equal  to  zero. 

Figure  4  illustrates  the  effect  of  mean  pressure  on  the  stress-strain 
curve  for  a3  =0,  a5  =  0.5,  and  n  =  0.5.  As  shown  in  the  figure,  the  strength 
of  the  material  increases  with  an  increase  in  mean  pressure. 

The  effect  of  the  strain  gradient  on  the  stress-strain  curve  is  shown  in 
Figure  5  for  aj  =  0,  a5  =  0.5  and  n  =  0.5.  One  of  the  assumptions  used  in  the 
current  model  is  displayed  in  this  figure,  i.e.,  the  stress  is  a  function  of 
both  strain  and  strain  gradient.  The  strength  of  the  material  is  reduced  with 
an  increase  in  strain  gradient. 

To  conclude  this  section.  Figure  6  shows  the  loading,  unloading,  and 
reloading  of  the  model  by  changing  the  sign  of  Ae12  for  al  =  0,  a3  =  0,  ' 
a5  =  0.5  and  n  =  0.5. 

The  above  figures  show  that  the  current  model  can  represent  many  response 
features  that  are  typical  for  geological  materials  (Refs.  1,  3,  and  7).  If 
larger  values  of  ae12  are  used,  the  response  becomes  less  smooth,  but  the 
essentia]  features  are  retained. 
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IV.  COMPARISONS  WITH  OTHER  NUMERICAL  DATA 

DYNAMIC  RESPONSE  OF  THE  INTERFACE  UNDER  SHEAR 

The  constitutive  relation  for  a  material  point  in  soil  has  been  estab¬ 
lished.  In  order  to  obtain  the  global  response  of  the  soil-concrete  interface 
under  dynamic  shear  loading,  the  finite  element  method  (Ref.  15)  together  with 
a  time  integrator  (Ref.  16)  can  be  used.  These  two  discretization  processes 
provide  approximate  solutions  of  complex  problems  for  which  analytical  solu¬ 
tions  might  not  be  available. 

In  this  section,  finite  element  approximations  over  a  given  geometrical 
domain  are  described  first.  Then,  the  field  variable  is  discretized  in  the 
time  domain  by  using  a  Newmark  integrator  in  time. 

For  the  soi 1 -concrete  interface  shown  in  Figure  1,  the  notations  of  engi¬ 
neering  shear  stress  (t)  and  strain  (y)  are  used  for  convenience.  It  is 

assumed  that  T  =  ai2,  and  that  the  displacement  (u)  is  a  function  of  (x)  only, 
and  in  a  unique  direction,  perpendicular  to  (x). 

With  these  assumptions,  it  follows  that 

Y  -  u,x  (24) 

and  the  one-dimensional  wave  equation  for  shear  is  given  as 

t,x  -  pu,u  (25) 

where  p  is  the  mass  density.  ! 

j 

i 

In  the  elastic  regime,  stress  is  related  to  the  displacement  throuah  the  \ 

; 

relation:  ! 


t  =  Gy  =  - 1 -  u  (26) 

2(1  +  v)  X 

Plasticity  theory  is  used  to  determine  stress  in  the  hardening  and  softening 
regimes. 


25 


For  the  soil  bar  of  length  (L),  the  variational  formulation  of  Eq.  25  is 


L 

/  (x,  -  p  u,  )  w  dx  =  0  for  all  w  (27) 

X  tc 


where  w  is  the  test  or  weighting  function.  After  performing  an  integration  by 
parts,  it  follows  that 


L  L 

/  w,x  t  dx  +  /  p  w  u,tt  dx 

o  o 


for  all  w 


(28) 


where  f.  and  f  are  the  external  loads  applied  on  the  ends  of  the  soil  bar, 
L  o 

x  =  L,  and  x  =  0,  respecti vely. 

Divide  the  soil  bar  into  elements  of  equal  length  (h)  as  shown  in  Fig¬ 
ure  7  and  use  linear  element  shape  functions  i.e.. 


4-1  *  1  -  n 


(29) 


2 


a  n 


(30) 


where 


n 


(31) 


Then  the  integrals  in  Eq.  28  can  be  converted  to  a  sum  with  the  result  that 

T  ( h  }  pe  we  u®  dn  +  /  w?  xe  dn')  =  —  w.  +  —  w  (32) 

0  o  n  /  A  L  A  0 

Q 

for  all  w,  where  N  is  the  number  of  elements  and  p  is  the  element  mass  den¬ 
sity.  The  global  functions  u,tt,  t,  w,  and  w,x  are  computed  element  by 
element  in  the  domain  of  the  problem  0  £  x  _<  L.  In  an  arbitrary  element  (e), 
these  functions  are  interpolated  as  follows: 


Figure  7.  Finite  element  mesh  of  the  soil  bar  and  the  shape 
function  used. 
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e  el  e2 

u‘tt  =  u*tt  +  u,tt 

(33) 

e  el  ,  ,e2, 

t  =  t  (pj  »  t  t|>2 

(34) 

e  el  e2 

w  1  w  1  fj  +  w  f  2 

(35) 

w>»  a  w*‘  »l,n  *  1,62  *2,(1 

(36) 

where  the  superscript  el  denotes  the  left  node  of  each  element,  and  e2,  the 
right  node. 


If  the  consistent  element  mass  matrix  is  defined  to  be 


Mu  * p  h  £  *i  *j 


dn 


and  the  element  internal  force  vector  is 


„Ie  \  e 

Ti  "  J  ♦i  n  1  dn 

1  o  i  »n 


then,  an  alternate  form  of  Eq,  32  is 


f  (  f  wei  M®  u  ej  +  I  f Ie  wei  }  =  fl  w  +  f°  w 

k  \i  ,j=i  ’u  i =1  1  J  VL  ~a  0 


for  all  w. 


(37) 


(38) 


(39) 


Now  relate  the  global  nodal  values,  w^,  tt>  and  the  element  nodal  or 
wel ,  utjt,  as 


local  values,  we1 ,  utjt»  as  follows: 


wk  =  W' 


=  w 


k,tt 


=  „(k- D2  .  „kl 


u; 


tt 


=  u’tt 


(40) 

(41) 


where  the  subscript  k  refers  to  the  node  number,  and  the  two  superscripts 
denote  the  element  number  and  the  element  node  number,  respecti vely.  Assume  p 
is  constant  for  each  element.  Then,  with  M  representing  the  global  mass 
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matrix  and  f . ,  the  global  force  vector,  the  variational  formulation  of  the 
equation  of  motion  can  be  written  in  indicial  notation  as 


Mij  Uj,tt  “  fi 


where 


1/3  1/6  0  0 

1/6  2/3  1/6  0 
01/60  0 
0  0  0. 

»  •  •  » 

0  0  0  . 


-  (t12  +  t21) 

-  (t22  +  T31) 


Equation  42  represents  a  set  of  ordinary  differential  equations  defined 
in  the  time  domain.  The  global  force  vector  is  the  sum  of  the  external  force 
and  internal  force  vectors,  through  which  the  constitutive  equation  sub-rou¬ 
tine  is  used.  To  solve  this  set  of  equations,  a  time  integrator  is  used.  To 
make  the  system  explicit  in  time,  a  diagonal  mass  matrix  is  required.  The 
diagonal  mass  matrix  is  constructed  by  summing  each  row  of  the  global  mass  . 
matrix  as  shown  in  Eq.  43,  and  putting  the  sum  of  each  row  into  the  corre¬ 
sponding  diagonal  position.  Although  the  time  step  must  be  chosen  according 
to  the  stability  condition,  no  matrix  storage  is  required,  and  calculations 
using  the  explicit  method  are  much  simpler  than  those  using  the  implicit 
method. 

Introduce  the  discrete  values  of  time, 

tn  =  ns  for  a  fixed  time  step  s  (45) 


>  VrVjTX  K>_  „>  '».■*  >; 


where  the  superscript  n  is  an  integer  denoting  the  discrete  time,  and  let  u^ 
vn,  and  a"  be  approximations  to  displacement  ui (tn),  velocity  vi (t  ),  and 
acceleration  u.  tt(tn)  at  node  i  respectively,  A  particular  form  of  tne 
Newmark  integrator  is 


n+1  n  ,  n  ,  s2  n 
ui  s  ui  +  s  vi  +  -  ai 


n+l  n  , 
vi  *  v.  +  s 


f\  ,  n  n+i 
(1  -a)  ai  +  a 


To  introduce  numerical  damping,  a  value  of  a  =  0.8  was  used.  In  addi¬ 
tion,  Eq.  42  is  approximated  with  the  equation 

M,,a"+1.f"+1  (48) 

Hence,  for  given  initial  and  boundary  conditions,  the  field  variable  u. 
at  each  node  can  be  found  without  determining  the  tangent  stiffness  matrix. 
With  the  use  of  the  constitutive  subroutine  mentioned  earlier,  the  algorithm 
for  obtaining  the  dynamic  response  of  the  soil  is  as  follows: 

Step  1  —  Input  material  constants  and  control  parameters. 

Step  2-- Input  initial  and  boundary  conditions. 

Step  3 — Start  time  integration  with  time  step  (s): 

.n  n* 1 

t  S  t  +  S 


Step  4 — Calculate  displacements  at  each  node  (i) 


n  n-1  .  n-i  s2  .n_i 

u  =u  +  s  v  +  —  a 

i  i  i  2  i 


i  =  1,  2 .  nn 


and  replace  Uj  at  the  end  of  the  soil  bar  by  the  given  boundary 
condition. 


Step  5--Find  total  strains  and  strain  increments  at  each  element  (i)  by 


using 


un  ,  -  u" 

Y"  =  111 .  1  i 

Yl  h 


1,  2,...,  nn  -  1 


,  n  n”1 
4Yf  -  tj 

Step  6— Compute  inelastic  strain  gradients  at  each  element  (i)  according 
to  the  following  relation: 

,  i  »n  ,  i  ,n 

i  n  -  ( Y -j ) 

(§  )"  *— i  =  1,  2 .  N-l 

h 

Step  7— Call  constitutive  subroutine  to  get  the  stresses  (t?). 

Step_8— Calculate  the  force  vector  at  each  node  (i)  by  means  of  Eq.  44. 
Step  9— Calculate  accelerations  at  each  node: 


a"  *  _L  f" 
1  “if  1 


where  is  the  diagonal  element  of  the  diagonal  mass  matrix. 


Step  10--Calcu1ate  the  velocity  vector: 


vj.vj-1  ♦.[(!-.)  a--1  *  -  «7] 


L — If  the  results  t(x,t)  and  Y(x,t)  are  needed— print  out,  then  go 


to  Step  3. 


To  verify  the  current  approach,  comparisons  with  other  existing  data  were 
required.  Existing  experimental  data,  however,  do  not  show  the  relation 


between  the  size  of  the  localization  zone  and  the  applied  load.  An  alternate 
theoretical  approach  (Ref.  12)  with  a  nonlocal  continuum,  called  an  imbricate 
continuum,  illustrates  numerically  the  dynamic  response  of  a  one-dimensional 
bar  subjected  to  tensile  waves.  Two  tensile  waves  are  initiated  at  the  ends 
of  the  bar,  then  the  waves  meet  together  at  the  middle  point.  That  point  will 
be  the  first  to  go  into  the  postpeak  regime.  Other  material  points  of  the  bar 
near  the  middle  point  will  also  get  into  the  postpeak  regime  because  of  the 
nonlocal  approach.  Thus,  the  strain  softening  accompanied  by  localization  can 
be  simulated.  In  the  imbricate  continuum  model,  a  bilinear  stress-strain 
relation  is  adopted.  Although  the  problem  of  the  soil -concrete  interface 
seems  to  be  different  from  that  in  Reference  12,  the  basic  features  are  the 
same  in  both  cases  if  the  middle  point  of  the  one-dimensional  bar  is  changed 
to  a  rigid  constraint,  and  only  half  the  bar  length  is  considered.  Therefore, 
the  numerical  solutions  of  Reference  12  can  be  used  to  corroborate  the  inter¬ 
face  model  by  taking  the  concrete  part  as  a  rigid  body  and  applying  the  cur¬ 
rent  constitutive  model  to  the  soil  part,  as  shown  in  Figure  7. 

By  dividing  the  soil  bar  of  length  (L  *  1)  into  20  elements,  using  the 
time  step  s  =  0.2  h  *typ/G,  and  assigning  to  the  other  parameters  the  follow¬ 
ing  values: 


e[  =  0.16 


e  =  10- 


=  nr 


at  =  0 


a2  =  0.1  a3  =  0.08 


a4  -  0.1  /T  a5  =  0.6 


n  =  0.4 


an  approximate  solution  is  obtained.  Results  are  shown  in  Figures  8  to  19. 

For  purely  elastic  behavior,  Figure  8  shows  the  propagation  of  the  stress 
wave  at  various  times  for 


6  *  -  =1  p  =  1 

2(1  +  v) 


fL  =  0.4  HCt] 


Total  stress 


X-coordinate 


Figure  9.  Spatial  distribution  of  stress  at  various  times 
(from  Ref.  12). 
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Displacement 


X-coordinate,  £■ 


Figure  14.  Spatial  distribution  of  strain  at  various  times  using 
the  current  approach. 


Displacement 


Shear  strain 


at  t  s  1.7  for  different  meshes. 


The  material  parameters  (G  and  p)  are  chosen  so  that  the  wave  speed  equals 
unity.  Except  for  numerical  dispersion  and  dissipation,  the  approximate  solu¬ 
tion  is  a  reasonable  representation  of  the  exact  shear  wave. 

Figures  9  to  14  show  the  inelastic  solutions  and  the  one-to-one  compari¬ 
sons  with  numerical  solutions  of  Reference  12  for 

G  * - l -  *  2  p  =  2  f,  =  0,9  HCtJ 

2(1  +  v) 

L  =  0.5  flf  L0  -  fT 
o  P 

From  these  figures,  it  is  apparent  that  the  tendencies  of  the  dynamic 
responses  of  the  two  problems  are  alike,  but  with  the  simpler  computational 
approach  given  by  the  current  model.  When  the  shear  wave,  which  propagates 
from  the  right  end  of  the  soil  bar,  reaches  the  interface,  a  reflection  occurs 
and  the  stress  increases  to  the  limit  state  which  initiates  softening.  As  the 
figures  show,  the  region  next  to  the  interface  is  characterized  by  large 
strain  and  low  stress,  which  is  believed  to  be  representati ve  of  what  occurs 
at  an  interface. 

Figures  15  and  16  show  the  stress-strain  responses  at  t  *  1.7  of  ele¬ 
ment  1  and  element  20,  which  represent  points  near  the  interface  and  near  the 
free  end  of  the  soil  bar.  Figures  17  to  19  show  the  convergence  of  the  numer¬ 
ical  solution  by  using  different  finite  element  meshes. 

Although  there  is  some  numerical  noise  in  the  numerical  solution,  the 
localization  zone,  i.e.,  the  large  deformation  band  with  stress  being  reduced, 
can  be  predicted  near  the  soil -concrete  interface.  In  addition,  it  appears 
that  good  results  can  also  be  obtained  by  using  a  coarse  mesh. 

STATIC  RESPONSE  OF  THE  INTERFACE  UNDER  SHEAR 

The  global  response  of  the  soi 1 -concrete  interface  under  static  shear 
loading  is  investigated  in  this  final  subsection.  The  problem  geometry  is 
shown  in  Figure  1.  The  soil  bar  is  divided  into  10  linear  elements,  and  the 
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load  is  a  monotonical ly  increasing  displacement  at  x  -  L.  The  shear  strength 
of  the  interface  is  assumed  to  be  slightly  smaller  than  other  parts  of  the 
soil  bar  so  that  the  interface  will  initiate  failure. 

The  algorithm  used  for  getting  the  static  response  of  the  interface  is 
outlined  as  follows: 

Step  1-- Input  material  constants,  boundary  conditions,  and  control 
parameters. 

Step  2--Monotonical ly  increase  the  strain  in  the  first  element,  and  find 
the  corresponding  increment  of  stress,  which  must  be  the  same  for  every 
element. 

Step  3--Ca11  the  constitutive  equation  subroutine  for  each  element.  If 
the  results  are  required,  print  them  out. 

Step  4--If  the  interface  is  still  in  the  prepeak  regime,  go  to  Step  2. 
Otherwise,  go  to  Step  5. 

Step  5—  If  the  previous  value  of  stress  for  any  element  is  larger  than 
the  current  peak  stress  corresponding  to  that  element,  go  into  the  softening 
regime.  Otherwise,  use  elastic  unloading. 

Step  6 — Go  to  Step  2  for  the  next  increment  of  strain  in  the  first 
element. 

Figures  20  to  27  show  the  numerical  solutions  for  the  following 
parameters : 


0.5  fT 

Lp = 

rr 

n  =  0.5 

- £ -  a  2 

2(1  +  v) 

-i 

eL 

=  0.16 

°LS  = 

10-** 

a2  - 

0.1 

a3  =  0.1 

o.i  rr 

Ay  = 

0.05 

a  1  =  0 
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Figure  20. 
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Figure  24.  Shear  stress-strain  response  of  element  6. 
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r  stress 


Figure  25.  Shear  stress-strain  response  of  element  10. 


52 


Shear  strain 
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Figure  20  shows  the  strain  distribution  of  static  response  in  the  elastic 
regime  where  the  strain  distribution  is  uniform.  Figure  21  shows  the  strain 
distribution  of  the  static  response  in  the  postpeak  regime.  Note  that  the 
strain  distribution  is  not  uniform,  and  that  the  localization  zone  occurs  near 
the  interface.  Figure  22  shows  the  stress  distribution  that  corresponds  to 
the  strain  distribution  of  Figure  21.  Since  there  is  no  distributed  applied 
force,  the  stress  distribution  is  always  uniform  under  static  loading.  Fig¬ 
ures  23  to  25  show  the  stress-strain  responses  of  element  1,  element  6,  and 
element  10,  respectively.  Element  1  represents  the  region  near  the  interface, 
and  element  10,  the  opposite  end  of  the  soil  bar.  After  the  peak  state,  ele¬ 
ment  6  unloads  elastically  at  first,  and  then  goes  into  the  softening  regime. 
Figure  26  shows  the  convergence  of  strain  distributions  for  different  finite 
element  meshes.  Figure  27  shows  the  effect  of  strain  gradients  on  the  static 
response  of  the  soil  bar  as  reflected  through  various  values  of  a3.  For  a3  = 
0,  i.e.,  local  continuum  approach,  the  localization  zone  cannot  be  predicted. 
This  is  why  a  classical  iocal  constitutive  model  cannot  simulate  the  soil- 
concrete  interface. 

Figure  28  shows  the  relations  between  shear  force  and  displacement  under 
the  effect  of  different  mean  pressures,  aLP,  for 

e^  =  0.4  as  =  0.5  a3  *  0.02 

with  other  parameters  being  unchanged.  According  to  experimental  results 
reported  by  Drumm  and  Desai  (Ref.  3),  basic  features  of  the  static  response 
near  the  interface  are  simulated  with  this  approach. 
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Displacement, 


Figure  28 
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Shear  stress  versus  the  displacement  of  the  end  of  the  soil 
bar  for  various  values  of  mean  pressure. 


V.  CONCLUSIONS 


It  was  assumed  that  no  relative  slip  occurs  at  the  contact-  surface 
between  soil  and  concrete.  Instead,  a  shear  band  appears  adjacent  to  the 
interface.  The  shear  band  is  just  the  physical  representation  of  softening 
and  localization.  The  interface  serves  as  an  initiator  of  strain  softening. 

Through  numerical  testing  and  indirect  comparisons,  it  was  concluded  that 
the  new  nonlocal  constitutive  model  can  simulate  the  softening  behavior  of 
different  materials  by  selecting  appropriate  values  of  material  parameters. 

The  approach  displays  the  same  features  as  those  in  Reference  12,  but  with  a 
much  simpler  computational  procedure.  Numerical  solutions  under  dynamic  and 
static  shear  loadings  illustrate  the  perceived  characteristics  of  a  soil- 
concrete  interface. 


Because  of  the  difficulties  in  measuring  the  localization  zone  with  high 
strain  gradients,  existing  experimental  data  are  not  satisfactory .  Once  good 
experimental  data  are  available,  it  will  be  necessary  to  verify  the  model,  not 
only  in  one-dimensional ,  but  also  in  two-  and  three-dimensional  cases  such  as 
curved  interfaces.  Further  progress  on  understanding  interface  phenomena  will 
require  a  combination  of  experimental  and  theoretical  investigations.  The 
research  reported  here  indicates  that  a  nonlocal  constitutive  model  is  a 
promising  theoretical  approach. 
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LIST  OF  SYMBOLS 


i 

L 


6EL 

E 


f. 


VfL 


F 


6 


H 

-i 

9 

h 


Total  strain  tensor 

Inelastic  deviatoric  strain  tensor 

Elastic  strain  tensor 

Inelastic  strain  tensor 

Inelastic  strain  invariant 

Limit  value  of  e1 

Shear  strain  at  the  elastic  limit 

Young's  modulus 

Global  force  vector 

External  loads 

Element  internal  force  vector 
Yield  surface  function 
Shear  modulus 

Strain  gradient  softening  function 
Strain  hardening  and  softening  function 

Inelastic  strain  gradient  invariant 
Length  of  a  finite  element 
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Length  of  a  soil  bar 


Elastic  limit  of  yield  surface 
Peak  value  of  yield  surface 
Element  mass  matrix 
Global  mass  matrix 
Material  parameter 
Mean  pressure 

Time  step  in  time  integration 
Global  displacement  vector 
Global  velocity  vector 
Weighting  function 
Engineering  shear  strain 
Total  strain  increment  tensor 
Elastic  strain  increment  tensor 
Inelastic  strain  increment  tensor 

Poisson's  ratio 
Mass  density 
Stress  tensor 
Engineering  shear  stress 
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